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Abstract 

In  the  discipline  of  non-Newtonian  materials,  the  ability  to  control 
viscoelastic  stresses  is  very  desirable  in  ascertaining  important  prop¬ 
erties  of  the  influenced  materials.  We  apply  the  nonlinear  geometric 
control  theory  to  examine  the  controllability  of  various  popular  consti¬ 
tutive  models  with  imposed  homogeneous  extensional  flow.  The  subse¬ 
quent  constitutive  laws  considered  here  include  the  Phan-Thien-Tanner 
model,  the  Johnson-Segalman  model,  the  Giesekus  model  and  the  Doi 
model.  This  paper  provides  the  first  analysis  on  the  effect  of  extensional 
flow  on  these  models. 


1  Introduction 

For  the  design  of  certain  materials,  controllability  of  viscoelastic  fluids  is 
a  signihcant  characteristic.  In  practice,  the  shape  of  an  extrudate  can  be 
controlled  by  varying  the  size  or  shape  of  an  orihce,  whereas  the  advance  of 
the  free  surface  can  be  controlled  by  varying  the  inflow  into  a  mold. 

Previous  theoretical  studies  from  Renardy  [[Renardy,  2005a],  [Renardy,  2005b], 
[Renardy,  2007]]  investigated  the  controllability  of  flows  of  linear  viscoelas¬ 
tic  fluids  for  the  multi-mode  Maxwell  models,  the  controllability  of  the  ho¬ 
mogeneous  shear  flow  of  viscoelastic  fluids  with  several  different  constitutive 
models,  and  the  controllability  of  nonhomogeneous  shear  flow  of  an  upper 
convected  Maxwell  fluid.  Very  recent  studies  have  examined  the  control¬ 
lability  of  the  upper  convected  Maxwell  model  under  various  homogeneous 
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flows  [[Zhou  et  ai,  2007]].  Furthermore,  the  reachable  set  for  the  upper  con- 
vected  Maxwell  model  under  imposed  extensional  flow  was  precisely  specihed 
in  [[Zhou  et  ai,  2007]]  even  though  it  is  usually  very  challenging  for  general 
cases.  In  most  of  these  studies  the  state  of  the  system  is  characterized  by  the 
viscoelastic  stresses  while  the  control  input  is  in  the  form  of  the  flow  rate  or 
the  body  force. 

The  primary  purpose  of  this  paper  is  to  further  these  earlier  studies  to 
consider  the  controllability  of  three  different  model  systems  under  the  imposed 
homogeneous  extensional  flow.  Like  shear  flow,  the  extensional  flow  is  also  very 
important  in  physical  applications.  For  example,  extensional  flow  can  be  used 
to  locally  approximate  the  flow  away  from  boundaries  in  extrusion  manufactur¬ 
ing  or  in  thin  him  and  sheet  manufacturing;  an  extension-dominated  how  oc¬ 
curs  in  industrial  wire-coating  processes  along  the  wire-coating  region  beyond 
the  die.  Following  earlier  studies,  we  restrict  our  attention  to  two-dimensional 
models.  This  is  mainly  due  to  its  mathematical  simplicity  and  is  also  physi¬ 
cally  motivated  by  monolayer  thin  hims.  Admittedly  two-dimensional  models 
are  ideal  situations  but  we  can’t  hope  to  understand  the  controllability  of  more 
realistic  cases  unless  we  understand  the  controllability  of  simpler  ones. 

For  reader’s  convenience,  we  now  present  on  overall  basic  dehnition  and 
description  of  weak  controllability  adopted  from  [[Isidori,  1995]]. 

Let 

m 

X  =  f{x) +  '^gi{x)ui  (1) 

i=l 

be  a  general  nonlinear  control  system,  where  x  is  the  state  variable,  and  Ui  G  M, 
i  =  1,  •  •  •  ,m,  are  the  control  variables.  Let  M  denote  the  manifold  of  state 
variables.  A  point  xi  in  M  is  reachable  from  a  point  xo  in  M  if  there  exist 
piecewise  continuous  input  functions,  Ui  =  ai(t),  snch  that  the  trajectory,  x(t), 
of  (1)  with  initial  state  Xq  reaches  Xi  in  hnite  time.  The  global  reachability  is 
nsnally  hard  to  prove  for  nonlinear  control  systems.  A  feasible  solution  is  to 
seek  weak  controllability.  The  system  (1)  is  weakly  controllable  within  some 
open  snbset  S'  C  M  if  for  each  point  xq  G  S',  there  is  an  open  neighborhood 
Uq  of  Xq  so  that  the  set  of  points  reachable  from  xo  along  trajectories  inside 
Uq  contains  at  least  an  open  snbset  of  M. 

In  (1),  f{x)  and  gi{x)  {i  =  are  vector  helds  in  M.  Under  the 

Lie  bracket  operation,  [f,g],  the  space  of  smooth  vector  fields  on  M  forms  a 
Lie  algebra.  This  Lie  algebra,  which  is  the  smallest  snbalgebra  containing  the 
vector  helds  /,  gi,  ■  ■  •,  gm,  is  called  the  Control  Lie  Algebra. 

A  powerfnl  snfhcient  condition  on  the  weak  controllability  of  a  nonlin¬ 
ear  control  system  is  the  contol lability  rank  condition  (CRC)  [[Isidori,  1995], 
[Hermann  and  Krener,  1977]].  Namely,  a  control  system  of  the  form  (1)  is 
weakly  controllable  on  an  open  set  S  if  it  satishes  the  controllability  rank 
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condition  on  S,  i.e., 

dim(Ac(x))  =  n  (2) 

for  all  X  E  S,  where  n  is  the  dimension  of  the  manifold  M,  and 

=  span{X(a;)|Xis  a  vector  field  in  the  control  Lie  algebra}. 

Note  that  the  big  advantage  of  the  CRC  is  that  it  is  purely  algebraic  and  it 
does  not  require  integrations  of  the  differential  equations. 


2  The  Phan-Thien- Tanner  Model 


The  Phan-Thien- Tanner  (PTT)  model  is  one  of  the  most  applied  differen¬ 
tial  type  of  the  nonlinear  viscoelastic  constitutive  equations.  It  contains  two 
parameters  that  control  nonlinearity  and  has  the  ability  to  £t,  to  some  extent, 
the  shear  and  elongational  properties  of  the  viscoelastic  materials  indepen¬ 
dently.  Thus,  the  Phan-Thien- Tanner  model  is  physically  more  realistic  than 
the  constant  shear  viscosity  Oldroyd-B  model  or  upper-convected  Maxwell 
model. 

The  Phan-Thien- Tanner  model  has  the  constitutive  equation 
[[Phan-Thien  and  Tanner,  1977]] 

t  -  (Vv)T  -  T(Vv)'^  +  AT  +  z/(trT)T  =  2/iD,  (3) 


where  T  is  the  stress  tensor,  v  is  the  velocity  vector,  Vv  is  the  velocity  gradient 
tensor,  A  is  the  relaxation  rate,  z/  is  a  constant,  the  notataion  “tr”  stands  for  the 
trace  of  the  tensor,  p  is  the  elastic  modulus  and  D  is  the  rate-of-deformation 
tensor. 

We  consider  2-D  homogeneous  viscoelastic  fluids  and  denote  the  stress  ten¬ 
sor  and  the  imposed  extensional  flow  with  rate  y{t)  by 


T 


Til  Ti2 
Ti2  T22 


v  =  (7(t)|,-7(t)|), 


(4) 


respectively,  where  Tn  is  the  first  normal  stress  difference,  T22  the  second 
normal  stress  difference  and  T12  the  shear  stress.  Suppose  the  control  input  is 
denoted  by  7  which  is  closely  related  to  the  velocity.  Then  the  general  dynamic 
problem  of  the  PTT  model  (3)  takes  the  form 

t  =  F(7(t),  T),  T(0)  =  To,  T(t;w)  =  Ti,  (5) 


where  Tq  and  Ti  are  the  given  initial  and  final  states.  The  state  of  the  system 
(5)  is  characterized  by  viscoelastic  stress  T  with  three  components  Tn,  T22 
and  Ti2. 
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For  the  extensional  flow  in  (4),  the  velocity  gradient  is 

'm  o' 

Vv=  ^  .  (6) 

2  . 

which  implies  that  the  rate-of-strain  tensor  is 

D  =  i|Vv  +  (Vv)T  =  M  J  ,  (7) 

Applying  the  velocity  gradient  (6)  and  the  rate-of-strain  tensor  (7)  to  the  PTT 
system  (3),  we  have  in  component  form  that 

^11  =  —[A  -|-  z^(Tii  -|-  T22)]Tii  +  {fi  +  Tii)'j{t) 

T22  =  —[A  -|-  viTll  +  T22)\T22  —  (h  +  T22)'i{t)  (8) 


Ti2  —  —[A  -|-  viTll  +  T22)\Ti2 

For  mathematical  convenience,  the  following  notation  is  introduced: 


X=  X2  =  T22  ■  (9) 

.  2:3  J  L  ^12  . 

Then  the  system  (8)  can  be  rewritten  as 

^  =  /(F)  +  ^(fK  (10) 

where 

L  +  xi 

f{x)  = -[\  +  u{xi+X2)]x,  g{x)  =  -{fi  +  X2)  ,u  =  y(t).  (11) 

_  0 

The  first  Lie  bracket  is  defined  by 

iyxi{xi  -  X2)  +  + 12x1  +  12x2) 

[f,9\='^9- f-^f-9=  iyx2ixi  -  X2)  -  +  uxi  +  12x2)  .  (12) 

12{XI-X2)X3 

Similarly,  the  second  Lie  bracket  is  given  by 

-iy\xi{xi  -  X2)  +  fi\[\  +  iy{xi  +  X2)] 

[/>  [/>  9]]  =  V[/,  g]-f-Vf-[f,  ^  =  -iyXx2{xi  -  X2)  -  /iA[A  +  u{xi  +  X2)]  . 

-iy\{xi  -  X2)X3 


(13) 
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Utilizing  the  vectors  from  g{x),  the  two  Lie  brackets  (12)  and  (13),  one  con- 
tructs  a  matrix.  By  investigating  the  rank  of  this  matrix,  the  weak  control¬ 
lability  can  be  characterized.  More  specihcally,  if  the  matrix  is  full  rank  or 
nonsingular,  then  the  system  is  weakly  controllable;  if  the  matrix  is  not  full 
rank  or  singular,  then  this  approach  offers  no  insight  to  the  weak  controllability 
of  the  system  and  one  has  to  seek  more  advanced  tools. 

We  now  compute  the  determinant  and  obtain 


det 


9,  [/,  [/,^]  =  2/iAz/a;3[A  +  +  X2)]ixi  -  X2f 


(14) 


The  determinant  does  not  equal  zero  when  x^  ^  0,  X1+X2  7^  — A/z/  and  xi  7^  X2. 
It  can  be  concluded  that  the  system  of  the  PTT  model  under  extensional  flow 
satishes  the  CRC  and  thus  is  weakly  controllable  when  T12  7^  0,  Tn  -|-  T22  7^ 
—  A/z/  and  Tn  7^  T22.  This  result  can  be  summarized  geometrically  as  follows. 
Given  =  {{xi,X2iX^)\xi,X2,Xz  G  i?},  there  exist  three  surfaces: 


51  =  {(a;i,a;2,a;3)|a;3  =  0} 

52  =  {(a;i,a;2,a;3)|A -f  z/(a;i  +  X2)  =  0} 
S'3  =  {(a;i,a;2,a;3)|a;i  =  X2} 


such  that  the  PTT  model  (3)  under  extensional  flow  is  weakly  controllable  at 
all  points  in  \  {S'!  U  5*2  U  S'3}. 


3  The  Johnson- Segalman  model 

The  Johnson-Segalman  model  characterizes  the  behavior  of  non-Newtonian 
fluids,  including  special  cases  of  Newtonian  and  Maxwell  fluids.  Additionally, 
it  is  a  viscoelastic  fluid  model  which  was  developed  to  allow  non-affine  defor¬ 
mations  [[Johnson  and  Segalman,  1977]]. 

The  Johnson-Segalman  (JS)  model  gives  a  viscoelastic  constitutive  equa¬ 
tion  in  the  form 

[(Vv)T  +  T(Vv)'^]  [(Vv)^T  +  T(Vv)]  +AT  =  2/iD.  (15) 

Here  a  is  a  parameter  describing  polymer  slip  where  —  1  <  a  <  1;  When  a  =  1, 
the  model  (15)  reduces  to  the  Oldrody-B  model. 

Analogous  to  the  previous  model,  we  now  apply  the  gradient  of  the  homo¬ 
geneous  extensional  flow  (6)  to  the  JS  system  (15)  which  leads  to 

Til  =  +  (h  + 

T22  =  -AT22  +  aT22)j{t)  (16) 


Ti2  —  — AT12 
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Note  that  the  third  equation  in  (16)  can  be  solved  exactly  such  that  Tu  = 
7"i2(0)e“^*.  Consequently,  the  Johnson-Segalman  model  (16)  is  not  weakly 
controllable  under  extensional  flow.  Nonetheless,  the  state  space  has  a  stable 
invariant  subspace  T12  =  0.  It  is  in  this  subspace  that  all  trajectories  of  the 
system,  under  any  control  input,  asymptotically  move  toward  the  subspace 
Ti2  =  0.  As  a  result,  the  decisive  behavior  of  the  control  system  (16)  can  be 
characterized  by  a  reduced  subsystem  on  this  stable  subspace: 


Til  —  “ATii  +  (/i  +  aTii)7(t), 

T22  =  — AT22  —  {fi  +  aT22)'j{t), 


(17) 


which  can  be  put  in  a  compact  form 


dx 

dt 


f{x)  +  g{x)u, 


where 


X  = 


xl 

'Til  ■ 

.  ^2  _ 

.  ^22 

-\x,  g{x) 


We  compute  the  Lie  bracket  and  And 


p  CtXi 

-(/i  +  0x2) 


u  =  7(«). 

(18) 


A  /i 
—A  n 


It  follows  immediately  that 


det 


aXg,{x2  —  Xi). 


(19) 


The  determinant  does  not  vanish  when  xi  7^  X2-  Therefore,  the  subsystem 
(17)  of  the  Johnson-Segalman  model  under  extensional  flow  satisfies  CRC  and 
is  weakly  controllable  when  Tn  7^  T22.  In  geometrically  terms,  the  system  (17) 
is  weakly  controllable  at  all  points  in  \  Si  where 

=  {(a;i,a;2)|a;i,a;2  e  R},  Si  =  {{xi,X2)\xi  =  X2}. 

In  [[Zhon  et  ai,  2007]]  we  have  characterized  the  set  of  reachable  states  for 
a  snbsystem  of  the  npper-convected  Maxwell  model  under  extensional  flow. 
We  have  fonnd  that  for  the  system 

Til  =  -(A  -  7(f))Tii  +  /i7(f) 


T22  =  -(A  -h  S(t))T22  -  p 7(f) 


(20) 
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the  reachable  set  is  i?i  +  Rq  where  Ri  and  Rq  are  given  by 


'Tiijtf) 


T22(tf) 


R,  =  {  (Tn(t/),T22(t/)) 


+  i)-i 


>  e-%|./(^  +  l)(Sfl  +  l)  -  1] 


'Tiijtf) 


T22{tf) 


+  1)  -1 


Ro  -  <  (Tu(tf),T22(tf))  =  e~^b[y/('^iho)  > 

^  and  Tii(O)  =  T22(0),  Tii(tf)=T22(tf)  ^ 

Note  that  if  we  redehne  the  two  parameters  /i  and  yit)  as  /ineto  =  jnewit)  = 
aioidit),  then  the  system  (17)  has  the  same  form  as  (20).  Therefore,  the 
reachable  set  for  the  snbsystem  of  the  Johnson-Segalman  model  (17)  nnder 
extensional  flow  is  i?i  +  Rq  where  /i  is  replaced  by  /i/a. 

4  The  Giesekus  Model 

Another  typical  nonlinear  viscoelastic  flnid  model  is  the  Giesekus  model 
[[Giesekus,  1982]].  The  constitutive  relation  for  the  Giesekus  model  is  given 
by 

t  -  (Vv)T  -  T(Vv)^  +  AT  +  uT^  =  2/iD.  (21) 

Applying  the  gradient  of  the  homogeneous  extensional  flow  (6)  to  (21),  we 
obtain 

Tn  =  -ATn  -  +  T/^)  +  (/i  +  Tn)7(t) 

T22  =  -AT22  -  u{T^2  +  Ti2)  -  (h  +  T22)i{t)  (22) 

Ti2  =  —[A  +  Z/(Tii  +  T22)]Ti2 

which  in  turn  can  be  put  in  the  compact  form 

(1  IT  -* 

-  =  f{x)+mu.  (23) 


r  1  r  T  1  r  +  ^3) 

Xi  ill 

X=  X2  =  T22  ,  fix)  =  -Xx2  -  yixl  +  xl) 

^3  Ti2  r\  ,  ™  m™ 


■[A  +  z/(a;i  +X2)]x^ 


/i  +  Xl 

gix)=  -(/i  +  X2)  ,  M  =  7(t). 

0 
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The  two  Lie  brackets  are  computed  as  follows: 


—  xl)  +  A/i  +  2h'fj,xi 

UX2{—X2  —  2fi)  —  \fi  +  uxl  , 
^{xi  -  X2)X3 


and 


(25) 


[/,  [fj]]  = 


As  a  result,  the  determinant  is  found  to  be 


a;|)(— A  +  2/iz/)  +  juX{X  +  2uxi) 

-  a;|)(A  —  2yuz/)  —  faX{X  +  210x2) 
iyx3{xi  -  X2){2iyfi  -  A) 

(26) 


det 


=  2/ira3(a;i -a;2)(A  -  z//i)^(A  -  2z//i).  (27) 


Clearly,  the  determinant  does  not  vanish  when  0:3  7^  0,  A  7^  z/jU,  A  7^  2z//i  and 
xi  ^  X2.  We  conclude  that  the  system  of  the  Giesekus  model  under  extensional 
flow  satishes  the  CRC  and  thus  is  weakly  controllable  when  T12  7^  0,  Tn  7^  T22 
and  the  parameters  satisfy  the  constraints  A  7^  z//i  or  2z/yU.  This  result  can  be 
summarized  geometrically  as  follows. 

Dehne  two  surfaces 


Si  =  {{XI,X2,X3)\X3  =  0}, 
*S'2  =  {(a;i,a;2,a;3)|a;i  =  X2}. 


Then  the  Giesekus  model  (21)  under  extensional  flow  becomes  weakly  control¬ 
lable  at  all  points  in  \  {S'!  U  S'2}  provided  that  the  parameters  satisfy  the 
constraints  A  7^  z//i  or  2z//i. 


5  The  Doi  Model 

The  Doi  model  for  rodlike  liquid  crystal  polymers  in  a  solvent  is  well- 
known  for  its  capability  to  describe  both  the  isotropic  and  nematic  phases 
and  phase  transition  between  them.  A  fundamental  element  of  the  model 
is  the  single  molecule  orientation  distribution  function.  Interactions  between 
molecules  are  represented  by  a  mean-held  potential.  The  rodlike  molecules  are 
also  subject  to  Brownian  force  due  to  the  fact  that  they  interact  with  other 
rodlike  molecules  and  with  the  how.  Generally,  the  model  is  a  microscopic 
Smoluchowski  equation  or  Fokker-Planck  type  equation  for  the  dynamics  of  the 
orientational  distribution  function  coupled  with  a  macroscopic  hydrodynamic 
equation  [[Doi  and  Edwards,  1986]].  The  Smoluchowski  equation  depicts  the 
convection,  rotation  and  dihusion  of  the  rodlike  molecules. 
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The  full  Doi  orientation  tensor  theory  is  developed  after  the  kinetic  Smolu- 
chowski  equation  is  projected  onto  a  second-moment  description  using  various 
closure  rules.  The  major  ingredient  in  this  tensor  theory  is  the  second-moment 
tensor  which  describes  the  orientational  distribution  of  the  ensemble  of  rodlike 
macromolecules.  The  orientation  tensor  is  traceless  and  symmetric.  The  phys¬ 
ical  and  practical  significance  of  the  orientation  tensor  is  that  it  is  the  basis 
for  micro-scale  light  scattering  measurements  of  primary  axes  (“directors”), 
degrees  of  molecular  alignment  (“birefringence”),  and  normal  and  shear  stress 
measurements.  The  study  of  two-dimensional  liquid  crystal  polymers  has  been 
physically  inspired  by  monolayer  films.  Thin  films  of  liquid  crystal  polymers 
are  used  as  alignment  layers  for  liquid  crystal  displays  because  of  their  sta¬ 
bility  and  nonlinear  optical  properties.  A  lot  of  theoretical  and  experimental 
studies  have  been  devoted  to  the  two-dimensional  Doi  model  (for  example,  see 
[[Lee  et  ai,  2006]]  and  references  therein). 

The  two-dimensional  Doi  model  is  given  by  [[Lee  et  ai,  2006]] 

Q  =  f2Q  —  Qf2  -I-  a[DQ  -|-  QD]  +  aD  —  2aD  ;  Q(Q  +  2^  ~  6DrF{Q),  (28) 
where  Q  is  the  orientation  tensor  such  that 


Q 


Qn  Qi2 
Q21  Q22 


Qn  —  ~Q22-,  Qi2  —  Q21; 


n  is  the  vorticity  tensor  such  that  fl  =  |[Vv  —  Vv^];  a  is  a  dimensionless 
parameter  which  depends  on  the  molecular  aspect  ratio;  D  is  the  rate-of-strain 
tensor;  is  the  rotary  diffusivity;  F{Q)  is  dehned  by 


F(Q)  =  (1  -  y)Q  -  NQ^  +  NQ  :  Q(Q  +  i) 

and  is  a  dimensionless  concentration  of  nematic  polymers. 

As  with  the  previous  models,  the  gradient  of  the  homogeneous  extensional 
flow  (6)  is  applied  to  the  Doi  system  (28),  yielding 


Qn  —  —QDr  [(1  —  y)Qii  +  2NQii{Ql^  Q12)]  +  ~  ‘^QiQiiQ^ 

(29) 

Qi2  =  —QDrQi2  [(1  —  y)  +  2N{Ql^  Q12)]  ~  2aQuQi2'yit). 

Using  the  nematic  relaxation  time  scale  the  flow  held  and  orientation  dy¬ 
namics  of  (28)  can  be  non-dimensionalized.  The  key  dimensionless  parameters 
are  then  the  Peclet  number  Pe{t)  =  y{t)/D^  (the  shear  rate  normalized  with 
respect  to  nematic  relaxation  rate)  and  the  dimensionless  concentration  pa¬ 
rameter  N .  Rescaling  time  as  t  =  the  nematodynamic  model  (28)  in  the 
dimensionless  form  becomes 

Qn  =  —QQn  [1  ~  T  +  ‘^^{Qh  +  Q12)]  +  “  2(5ii)T’e(t), 

Qi2  =  —QQ12  [1  ~  Y  +  2A^((5ii  -|-  Q12)]  ~  2aQnQi2PG{t)- 


(30) 
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The  system  (30)  can  be  written  as 

^  =  fix)  +  gix)u, 

where 

^  ^  i  rM  ^  i  =  -6  [1  -  f  +  2iV(a^?  +  xD]  X, 

_  X2  j  I  V12  J 
a(|  —  2x1) 

(fix)  =  :  u  =  Pe(t). 

—2axiX2 

The  Lie  bracket  is  calculated  to  be 

■  6a  {2a;2[l  -  f  -  2N{xl  +  xj)]  +  |  -  f  +  iV(3a;2  +  xD]  ' 

[f,  ^  = 

12aa;ia;2[l  +  y  —  2N{x\  +  x\)] 

Finally,  we  have 

det  g,[f,g\  =l2afxiX2.  (32) 

The  determinant  does  not  equal  zero  when  xi  7^  0  and  X2  7^  0.  So  the  Doi 
model  (30)  under  extensional  flow  satisfles  CRC  and  is  weakly  controllable 
when  Qii  7^  0  and  Q12  7^  0.  In  terms  of  geometric  words,  the  system  (30)  is 
weakly  controllable  at  all  points  in  \  {Pi  U  P2}  where 

=  {(a;i,a;2)|a;i,a;2  G  R},  Pi  =  {(xi, a;2)|a;i  =  0},  P2  =  {(xi, a;2)|a;2  =  0}. 

6  Concluding  Remarks 

We  have  applied  the  controllability  rank  condition  to  the  vector  fields  in  the 
Phan-Thien-Tanner  model,  the  Johnson-Segalman  model,  the  Giesekus  model 
and  the  Doi  model  to  study  the  controllability  of  non-Newtonian  fluids  driven 
by  imposed  homogeneous  extensional  flow.  In  our  control  system,  the  state 
variable  is  the  stress  for  the  Phan-Thien-Tanner  model,  the  Johnson-Segalman 
model  and  the  Giesekus  model  and  the  orientation  tensor  for  the  Doi  model. 
The  extensional  flow  rate  corresponds  to  the  available  control.  We  have  derived 
sufficient  conditions  for  the  weak  controllability  of  each  model. 
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